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DISCLAIMER 



Although each program has been tested by its contributor, no 
warranty, express or implied, is made by the contributor or 
COMMON, as to the accuracy and functioning of the program 
and related program material, nor shall the fact of distribution 
constitute any such warranty, and no responsibility is assumed 
by the contributor or COMMON, in connection therewith. 



COMMON USERS GROUP PROGRAM REVIEW AND EVALUATION 
(fill out in typewriter, ink or pencil) 

Program No. Date 

Program Name: 

1. Does the abstract adequately describe what the program is and what Yes No_ 

it does ? 

Comment 

2. Does the program do what the abstract says? Yes No_ 

Comment 

3. Is the description clear, understandable, and adequate? Yes No_ 

Comment 

4. Are the Operating Instructions understandable and in sufficient detail? Yes No_ 

Comment 

Are the Sense Switch options adequately described (if applicable) ? Yes No_ 

Are the mnemonic labels identified or sufficiently understandable ? Yes No_ 

Comment 



5. Does the source program compile satisfactorily (if applicable) ? Yes Nq_ 

C om m ent 

jp^6. Does the object program run satisfactorily? Yes No_ 

%J Comment 

7. Number of test cases run . Are any restrictions as to data, 

size, range, etc. covered adequately in description? Yes No_ 

Comment 

8. Does the Program meet the minimal standards of COMMON? Yes No_ 

Comment 

9. Were all necessary parts of the program received? Yes No_ 

Comment 

10. Please list on the back any suggestions to improve the usefulness of the program. 
These will be passed onto the author for his consideration. 

Please return to: Your Name 

. Company 



Mr. Richard L. Pratt 

Data Corporation Address 

7500 Old Xenia Pike 

Dayton, Ohio 45432 Users Group Code 



|I|TH]S REVIEW FORM IS PART OF THE COMMON ORGANIZATION'S PROGRAM REVIEW AND 
^EVALUATION PROCEDURE. NONME MBERS ARE CORDIALLY INVITED TO PARTICIPATE 
IN THIS EVALUATION. 
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Explanation of Program 
The program fits a polynomial of order 9 to a set of data points, 
using a trucated series of Legendre polynomials. While the running time is 
somewhat longer than a least squares fit, the Legendre series seems to follow 
any "bumps" in the data curve more closely, and several people have preferred 
this effect to the "smoothing out of the bumps" which the least squares fit 
tends to do. 

One of the data input quantities needs a special explanation. A 

function f (x) can be represented by a Legendre series 
oo 2k+ 1 l 

f (x) kS> C k p k (x)# c k = — i f (x) p k (x) where p k w are 

Legendre polynomials. The higher order polynomials oscillate badly in the 
interval (-1,1), so the interval length chosen for the integral for C k must be 
small, in many cases many times smaller than the intervals between the given 
data points. Thus a quantity IJ is introduced to represent the number of sub- 
intervals between adjoining data points. Then f (x) is interpolated linearly at 
the ends of each subinterval. 

To take an example there should be a total of around 200 subintervals 
in the interval (- 1,1) for a reasonably good integration of P 9 (x) f(x). If there 
are, say, 12 data points, IJ should be about 20. 

Modifications or revisions to this program, as they occur 
will be announced in the appropriate Catalog of Programs' 
for IBM Data Processing Systems. When such an announce- 
ment occurs, users should order a complete new program 
from the Program Information Department, 







DATA INPUT FORMAT 
Punched Cards. Order of Data: N (Number of Data Points) IJ (Number of 
Subintervals Between Adjacent Data Points). 

Xj, F(X,), X 2 , F(X 2 ), X n , F(X r ), as closely as one can pack them on a 

card or cards. 

DATA OUTPUT FORMAT 

Punched Cards. 

If SSW1 is on, the Legendre series coefficients are punched. 

10 

If SSW2 is on, the coefficients of the power series f(x) 2= P w k -1 are 

k=l k 

punched. 

If Sense Switch 3 is on, the x values are "normalized" into the interval (-1, l), 
and these quantities are punched: x, w, s, d, for each data point, 
where, x = abscissa of data point 

w = normalized value of x 

s = value of truncated Legendre series 

d = difference between s and f, the ordinate of data point 
If Sense Switch 3 is off, the Legendre series in w is converted to a power series 
in x, and these quantities are punched: 

10 , . 

the coefficients r k of the power series f(x) = 2 r, x 

k= 1 k 

x, s,d for each data point, where x, s d are as defined above. 




MACHINE CONFIGURATION 
IBM 1620 with card reader- punch 60 K memory. 



PROGRAM INSTRUCTIONS 
Beginning of program: 07500 
Load card reader, push start. 

If Legendre series coefficients are desired, put Sense Switch 1 on. If 
normalized power series coefficients are desired, put Sense Switch 2 on. 
If summation of normalized power series is desired, put Sense Switch 3 on. 
This maybe necessary, because the regular power series is liable to overflow 
i? abscissas are large. 

When program is finished, new data is called for. 

Program running time in seconds: T « 1.5 x (Number of Data Points) x 
(Number of Subintervals between data points) = 1.5 x (N) x (IJ) seconds. 
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DO 



; INTEGRATION IS 
j PERFORMED ! 
I C k -2 .AAA.F.fyVl 

ri I$TEJ 



^ INTERVAL 
yiOT FI NISHED 



DO 59 



i k - ! 2 . ,10i INTERVAL FINISHED 

PROGRAM CONTINUES AT (2) <f . 

INTEGRATION FINISHED 
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DO 10 
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POWER SERIES 




POWER SERIES 
16,92 COEFFICIENTS ARE 
70,71 COMPUTED R. 

L.PQR ORIGINAL DATA._ 



ssw 3y> 



JQFF. " f: 
PUNCH 

\ .. . 

| SERIES IS SUMMED 
DO 80 jFOR EACH X IN DATA 
S 



... ... ON 

(IF OVERFLOW 
HAS OCCURRED) 



L 



/ssw 3 



MORE DATA £l 
RETURN TO BEGINNING 



__O.FF. 

PUNCH 
X,S,DIFF. 
FOR EACH X 



ON > 

(IF OVERFLOW ' 
HAS OCCURRED) 



END OF PROGRAM 
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SERIES IN W IS 
SUMMED FOR 
EACH X IN DATA S 



MORE DATA Si 
RETURN TO BEGINNING 



PUNCH 
X.W,S,DIFF. 

for EACH X 



END OF PROGRAM 
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POLYNOMIAL CURVE FITTING USING LEGENDRE POLYNOMIALS 



ct S™^'? 1 ! X(100) ' F(100 )' C ^°)' W ^00),V(9).P(10),A(9).B(9),R(10) 
51 KtAD,N, I J 

ET=IJ 

ET=1 .OOOOOOO/ET 
IJ-IJ+1 
DO 1 1=1 ,N 
READ,X(I),F(I) 
1 CONTINUE 
DO 20 J=1 , 10 
C( J) «.00000000 
20 CONTINUE 
DO 2 l**1,N 

, Sflii -2.0000000*( X ( I ) -X ( 1 ) ) / ( X ( N ) -X ( 1 ) ) -1 . 0000000 
Z CONT I NUE 
M=N-1 

DO 59 1=1, M 
Q«. 00000000 
T«ET*(W(I+1)-W(I)) 
TT«ET*(F(I+1)-F(I)) 
DO 60 J-1JJ 
V(1)«W( I )+T*Q 
FF-F(I)+TT*Q 
Q»Q+1. 0000000 
DO 3 K=2,9 
V(K)-V(K-1)*V(-1) 
3 CONTINUE 

P(1)*1 .0000000 
P(2)=V(1) 

P(3)=1 .5000000*V(2)-. 50000000 

P(k) -2 . 5000000*V ( 3 ) - 1 . 5000000*V ( 1 ) 

P( 5 )= t f.3750000*V(lf)-3.7500000*V( 2)+. 37500000 

I ) S < "1; 8 7 5 °000*V ( 5 j -8 . 75 00000*V ( 3 ) +1 . 8750000*V ( 1 ) 

o Z =1 5*S 37500 V 6 -19.687500*v(ij+6.5625000*V 2 -.31250000 

P(8)»26.812500*V(7H3.312500*V(5)+19.687500*V(3)-2^ 875000*V(1) 

^(9j:50-273^38*V(8)-93 .8l»375W«)^.1liM25^KM:ffl?5^(2) 

, IF (J-2) 4,5.5 
h AAA=.50000000*T 
GO TO 7 

5 AAA=T 

6 IF (J-IJ)7,8,8 
8 AAA=.50000000*T 

7 DO 9 K-1,10 
C(K)=C(K)+AAA*FF*P(K) 



9 CONTINUE 
60 CONTINUE 
59 CONTINUE 

DO 10 K=1,10 

E=K-1 

C(K)=.50000000*C(K)*(2.0000000*E+1 .0000000) 
10 CONTINUE 

I F( SENSE SWITCH 1)12,13 

12 PUNCH, C(1),C(2),C(3),C(lf) 
PUNCH, C( 5) ,C(6).C(7),C(8) 
PUNCH, C(9),C( 10) 

13 P(10)=9 4 ».960938*C(10) 

P 9)=50.273^38*C(9) n , ^ 

P(8)=26.812500*C(8)-201 .09380*C( 10) 

P(7)=1^^37500*C(7)-93.8^3750*C(9) 

P(6)=7.87500OO*C(6)-'t3.3l250O*C(8) Ht+0.76562*C(10) 

P(5)=*f. 3750000*C(5> -19. 687500*C( 7 )+5*t. 1<+0625*C(9) 

P(if)=2.5000000*C(^)-8.7500000*C(6)+19.687500*C(8)-36.093750*C( to) 

P(3)=1.5OOOOO0*C(3)-3.750O0O0'>C(5)+6.5625000*C(7)-9.8lf37500*C(9) 

P(2)=C(2)-1.5O00OO0*C(4)+t.8750000 ;t C(6)-2.l875000*C(8) 

P( 2 ) =P( 2)+2 .*f609375*C( 10) 

P ( 1 ) -C ( 1 > - . 50000000*C ( 3 )-»• ■ 3 75 00000*C ( 5 > - . 3 1 2 5 0000*C ( 7 ) 

P(1)=P(1)+.273 1 *3750*C(9> 

I F( SENSE SWITCH 2)15,16 

15 PUNCH, P(t).P(2),P(3),P( 1 ») 
PUNCH, P(5) ,P(6) ,P(7) »P(8.) 
PUNCH, P(9),P(10) 

16 AA«2. 0000000/ (X(N)-Xd>) 
BB=-AA*X( 1 )-1 .0000000 
A(1)=AA 

DO 70 1=2,9 
A(l)=AA*A(l-i) 

70 CONTINUE 
B(t)=BB 

DO 71 l=2 f 9 
B( l)=BB*B( 

71 CONTINUE 

I F( SENSE SWITCH 3)90,92 
92 R(1)=P(1)+P(2)^B(1)+P(3) V; B(2)!-P(4)^B(3)+P(5)*B(U)+P(6)*B(5) 
R(1)=R( 1)+P(7)*B(6)+P(8j*B-(7)+P(9)*B(8)+P(10)*B(q) 
R( 2 -P( 2 ) +? . 0000000*P ( 3 > * B( 1 ) +3 . 0000000*P( h ) *B ( 2 ) 
R ( 2 ) -R< 2 ) +k . 0000000* P ( 5 ) *B ( 3 ) +5 . 0000000*P ( 6 ) *B ( If ) 
R(2)=R(2)+6.0000000*P(7)» V B(5) + 7.0000000*P(8)*B(6) 
R( 2 ) =A( T ) *( Pv( ) +8 • 0000000*P( 9 ) *'B 1 7 ) +9 . 0000000*P( 1 0) *B ( 8 ) ) 
R( 3 ) =P ( 3 )+3 . 000000<)*P (k ) *B ( 1 ) +6 . 0000000*P( 5 ) *B ( 2 ) 
R( 3 ) »R ( 3 ) + 1 0. 000000* P( 6 ) *B ( 3 ) + 1 5 . 000000*P ( 7 ) *B ( «0 



7 
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R(3)=R(3)+21 .000000*P(8)*B(5)+28.000000*P(9)*B(6) 
R(3)=A(2)*(R(3)+36.000000*P(10)*B(7)) 
R(l+)=P(4)+^.0000000*P(5)*B(1 ) + 10.000000*P(6)*B(2) 
R(M=R(U)+20.000000*P(7)*B(3)+35.000000*P(8)*B(4) 

R(*O=A(3)*(R(2O+56.000000*P(9)*B(5)+8l*.000000*P(l0)*B(6)) 

R( 5 )=P( 5 ) +5 .0000000*P( 6 )*B( 1 )+1 5 . 000000*P( 7 )*B( 2 ) 

R(5)=R(5)+35.00O0OO*P(8)*B(3)+70.000000*P(9)*B(i») 

R(5)=A(U)*(R(5)+126.00000*P(10)*B(5)) 

R ( 6 ) =P ( 6 ) +6 . 0000000*P ( 7 ) *B ( 1 ) +2 1 . 000000*P ( 8 ) *B ( 2 ) 

R(6)=A(5)*(R(6)+56.O0000O*P(9)*B(3)+126.0OO0O*P(1O)*B(if)) 

R(7)=P(7)+7.0O0O00O*P(8)*B(l)+28.0O0OOO*P(9)*B(2) 

R(7)*A(6)*(R(7)+84.000000*P(10)*B(3)) 

R( 8 ) =A( 7 ) * ( P( 8 ) +8 . 0000000*P ( 9 ) *B( 1 ) +36 . 000000*P ( 1 0) *B ( 2 ) ) 
R(9)=A(8)*(P(9)+9.0000000*P(10)*B( 1)) 
R(10)-A(9)*P(tO) 
PUNCH, R( 1 ) ,R( 2) ,R(3) ,R(^) 
PUNCH, R( 5) ,R(6) ,R(7) ,R(8) 
PUNCH, R(9),R(10) 
I F( SENSE SWITCH 3)90,93 
93 DO 80 1=1 , N 
V(1)-X(I) 
DO 81 J=2,9 
V(J)-V(J-1)*V(1) 
IF(SENSE SWITCH 3)90,81 

81 CONTINUE 
S-R(|) 

DO 82 J=1.9 

S=S+R(J+1)*V(J) 

I F( SENSE SWITCH 3)90,82 

82 CONTINUE 
D=S-F( I ) 
PUNCH, X( I ) , S , D 

80 CONTINUE 

GO TO 51 
90 DO 17 1=1, N 

V(t)-W(l) 

DO 18 J=2,9 

V(J)-V(J-1)*V(1) 

18 CONTINUE 
S=P(1) 

DO 19 J=1,9 
S»S+P(J+1)*V(J) 

19 CONTINUE 
D=S-F( I ) 

PUNCH, X( I ) ,W( I ) , S,D 
17 CONTINUE 
GO TO 51 
END 
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